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Abstract 

The differential cross section da/dt for elastic p° photoproduction, jp — ► p°p (p° — ► 
7r + 7r~), has been measured in ep interactions at HERA. The squared four-momentum 
exchanged at the proton vertex, t, has been determined directly by measuring the mo- 
mentum of the scattered proton using the ZEUS Leading Proton Spectrometer (LPS), 
a large scale system of silicon micro-strip detectors operating close to the HERA pro- 
ton beam. The LPS allows the measurement of the momentum of high energy protons 
scattered at small angles with accuracies of 0.4% for the longitudinal momentum and 
5 MeV for the transverse momentum. Photoproduction of p° mesons has been inves- 
tigated in the interval 0.073 < \t\ < 0.40 GeV 2 , for photon virtualities Q 2 < 1 GeV 2 
and photon-proton centre-of-mass energies W between 50 and 100 GeV. In the mea- 
sured range, the t distribution exhibits an exponential shape with a slope parameter 
b = 9.8±0.8 (stat.)rtl.l (syst.) GeV -2 . The use of the LPS eliminates the contamination 
from events with diffractive dissociation of the proton into low mass states. 
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1 Introduction 



The physics of diffractive scattering processes has emerged as one of the most interesting topics 
of study in the early running of HERA. Up to now, the cross section for events in which 
the virtual photon diffractively dissociates into a vector meson or a generic state X has been 
measured in the HI and ZEUS experiments either by requiring a "large rapidity gap" between 
the proton beam direction and the most forward energy deposit recorded in the detector or by 
subtracting the non-diffractive background in a statistical way (lJ-0. Here we present the first 
cross section measurement at HERA in which diffraction is tagged by the detection of a high 
energy scattered proton, thereby eliminating contamination by events with dissociation of the 
proton. 

The measurement of the proton was performed using the ZEUS Leading Proton Spectrometer 
(LPS), which detects protons scattered at very small angles (^1 mrad). In this spectrometer, 
silicon micro-strip detectors are used in conjunction with the proton beam line elements to 
measure the momentum of the scattered proton. The detectors are positioned as close as the 
10<r envelope of the circulating proton beam (typically a few mm) by using the "Roman pot" 
technique ||. In the configuration used to collect the data presented here, the LPS consisted 
of a total of about 22,000 channels. 

This paper concentrates on the exclusive process 7p — > p°p in ep interactions at small photon 
virtualities (Q 2 ~ 0, the "photoproduction" region). This reaction is often called "elastic", in 
reference to the vector meson dominance model (VDM). Elastic photoproduction of p° mesons 
has been investigated in fixed target experiments at photon-proton centre-of-mass energies 
W < 20 GeV §-Q as well as at HERA energies, W « 100-200 GeV g |. The process has 



the characteristic features of soft diffractive interactions: the dependence of the cross section 
on W is weak, the dependence on t is approximately exponential, and the vector meson is 
observed to retain the helicity of the photon (s-channel helicity conservation). Here t is the 
squared four-momentum exchanged at the proton vertex. The data presented in this paper 
cover the kinematic range 50 < W < 100 GeV, Q 2 < 1 GeV 2 and 0.073 < |t| < 0.40 GeV 2 . 
Elastic events were selected by requiring that the scattered proton carry more than 98% of 
the incoming proton beam energy. The scattered positron was not detected; however, Q 2 was 
estimated using transverse momentum balance. 



2 Experimental set-up 
2.1 HERA 

The data discussed here were collected in 1994 at HERA which operated with 820 GeV protons 
and 27.5 GeV positrons (indicated in the following with the symbol e). The proton and positron 
beams each contained 153 colliding bunches, together with 17 additional unpaired proton and 
15 unpaired positron bunches. These additional bunches were used for background studies. 
The time between bunch crossings was 96 ns. The typical instantaneous luminosity was 1.5 x 
10 30 cm _2 s~ 1 and the integrated luminosity for this study is 898 ± 14 nb _1 . 
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2.2 The ZEUS detector 



A detailed description of the ZEUS detector can be found elsewhere JT2J. A brief outline of the 



components in the central ZEUS detector |13| which are most relevant for this analysis is given 
below, followed by a description of the Leading Proton Spectrometer. 



2.2.1 Central components and luminosity measurement 

Charged particles are tracked by the inner tracking detectors which operate in a magnetic field 
of 1.43 T provided by a thin superconducting coil. Immediately surrounding the beam pipe 
is the vertex detector (VXD), a drift chamber which consists of 120 radial cells, each with 12 
sense wires []I4" |. It is surrounded by the central tracking detector (CTD), which consists of 72 



cylindrical drift chamber layers, organised into 9 superlayers covering the polar angle region 
15° < 9 < 164°0 0. 

The high resolution uranium-scintillator calorimeter (CAL) |JT6[ consists of three parts: the 
forward (FCAL), the rear (RCAL) and the barrel calorimeter (BCAL). Each part is subdivided 
transversely into towers and longitudinally into one electromagnetic section (EMC) and one 
(in RCAL) or two (in BCAL and FCAL) hadronic sections (HAC). A section of a tower is 
called a cell; each cell is viewed by two photomultiplier tubes. The CAL energy resolution, as 
measured under test beam conditions, is ue/E = 0.18/v^E for electrons and &e/ E = 0.35/-\/E 
for hadrons (E in GeV). 

The Veto Wall, the C5 counter and the small angle rear tracking detector (SRTD) all consist 
of scintillation counters and are located at Z = —730 cm, Z = —315 cm and Z = —150 cm, 
respectively. Particles which are generated by proton beam-gas interactions upstream of the 
nominal ep interaction point hit the RCAL, the Veto Wall, the SRTD and C5 at different times 
than particles originating from the nominal ep interaction point. Proton beam-gas events are 
thus rejected by timing measurements in these detectors. 

The luminosity is determined from the rate of the Bethe-Heitler process, ep — > e'jp, where 
the photon is measured with a calorimeter (LUMI) located in the HERA tunnel downstream 



of the interaction point in the direction of the outgoing positrons [17 



2.2.2 The Leading Proton Spectrometer 

The Leading Proton Spectrometer |12| (LPS) detects charged particles scattered at small angles 
and carrying a substantial fraction, xl, of the incoming proton momentum; these particles 
remain in the beam pipe and their trajectory is measured by a system of position sensitive 
silicon micro-strip detectors very close to the proton beam. The track deflection induced by 
the magnets in the proton beam line is used for the momentum analysis of the scattered proton. 

The layout of the LPS is shown in Fig. |l]; it consists of six detector stations, SI to S6, placed 
along the beam line in the direction of the outgoing protons, at Z = 23.8 m, 40.3 m, 44.5 m, 
63.0 m, 81.2 m and 90.0 m from the interaction point, respectively. 

lr rhe coordinate system used in this paper has the Z axis pointing in the proton beam direction, hereafter 
referred to as "forward" , the X axis pointing horizontally towards the centre of HERA and the Y axis pointing 
upwards. The polar angle 9 is defined with respect to the Z direction. 
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Each of the stations SI, S2 and S3 is equipped with an assembly of six planes of silicon 
micro-strip detectors parallel to each other and mounted on a mobile arm, which allows them 
to be positioned near the proton beam. Stations S4, S5 and S6 each consist of two halves, each 
half containing an assembly of six planes similar to those of SI, S2, S3, also mounted on mobile 
arms, as shown in Fig. |^. Each assembly has two planes with strips parallel to the direction of 
motion of the arm, two planes with strips at +45° and two at —45° with respect to it; this makes 
it possible to measure the particle trajectory in three different projections in each assembly. The 
dimensions of the detector planes vary from station to station but are approximately 4x6 cm 2 . 
The pitch is 115 /im for the planes with vertical or horizontal strips and 115/v2 = 81 /im for 
the planes with ±45° strips. The distance along Z between neighbouring planes in an assembly 
is ~ 7 mm. The detector planes are mounted in each assembly with a precision of about 30 /im. 

The detector planes are inserted in the beam pipe by means of re-entrant "Roman pots" 
which allow the planes to operate at atmospheric pressure. A pot consists of a stainless steel 
cylinder with an open end away from the beam; the other end is closed. The silicon detector 
planes are inserted from the open end and are moved in until they are at about 0.5 mm from 
the closed end. The whole cylinder can be inserted transversely into the beam pipe. Figure [| 
illustrates the principle of operation. The walls of the pots are 3 mm thick, except in front of 
and behind the detector planes, where they are 400 /im thick; the thickness of the pot bottom 
walls facing the beam is also 400 /im. The vacuum seal to the proton beam pipe is provided by 
steel bellows. The pots and the detector planes are positioned by remotely controlled motors 
and are retracted during the filling operations of the collider to increase the aperture of the 
vacuum chamber; this also minimises the radiation damage to the detectors and the front- 
end electronics. In stations SI, S2, S3 the detector planes are inserted into the beam pipe 
horizontally from the outside of the HERA ring towards the centre. In stations S4, S5, S6, the 
detector planes in the two halves of each station independently approach the beam from above 
and from below. In the operating position the upper and lower halves partially overlap (cf. 
Fig. 0). The offset along the beam direction between the centres of the upper and lower pots 
is ~ 10 cm. Stations S5 and S6 were used in an earlier experiment at CERN and were adapted 
to the HERA beam line 



Each detector plane has an elliptical cutout which follows the profile of the 10<r envelope 
of the beam, where a is the standard deviation of the spatial distribution of the beam in the 
transverse plane. Since the 10<r profile differs from station to station, the shape of the cutout 
varies from station to station; in data taking conditions the distance of each detector from the 
beam centre is also different and ranges from 3 to 20 mm. Small variations of the detector 
positions from fill to fill are necessary during operation in order to follow the changes of the 
beam position and adapt to the background conditions. 

The detector planes are read out by two types of VLSI chips mounted on the detector support: 



a bipolar amplifier-comparator |19j followed by a radiation hard CMOS digital pipeline |20 
which operates with a clock frequency of 10.4 MHz, synchronous with the HERA bunch crossing. 
Each chip has 64 channels reading out 64 adjacent strips. The chips are radiation hard up to 
doses of several Mrad. 



A simplified diagram of the spectrometer optics is shown in Fig. |3], in which the beam line 
elements have been combined to show the main optical functions. Together with the HERA 
proton beam magnets, the six LPS stations form two spectrometers: 
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1. Stations SI, S2, S3 use the combined horizontal bending power of a septum magnet and 
three magnetic septum half-quadrupoles. SI, S2, S3 were not operational in 1994 and are 
not discussed further here. 

2. Stations S4, S5, S6 exploit in addition the vertical bending provided by three main dipole 
magnets (BU). These stations were used for the present measurement. 

The insertion of the detectors into the operating positions typically begins as soon as the 
beams are brought into collision. Among the conditions required prior to beginning the insertion 
are the following: (i) proton beam position as measured with the HERA beam position monitor 
next to S4 within 1 mm of the nominal position; (ii) background levels as measured in counters 
downstream of the main proton beam collimators, in the C5 counter and in the trigger counters 
of the Forward Neutron Calorimeter PI] (located downstream of S6 at Z m 109 m) stable 
and below given thresholds. About fifty minutes were necessary in 1994 to insert the detector 
planes. This and the fact that the beam conditions did not always allow safe insertion of the 
detectors results in the reduced value of the integrated luminosity available for this analysis 
with respect to other analyses of the ZEUS 1994 data. 

The strip occupancy during data taking, i.e. the average number of strips firing per trigger 
divided by the total number of strips, depended on the beam conditions but was typically less 
than 0.1%, with small contributions from noise and synchrotron radiation. 

The fraction of noisy and malfunctioning channels in 1994 was less than 2%; they were due 
to bad detector strips and dead or noisy front-end channels. The efficiency of the detector 
planes, after excluding these channels, was better than 99.8%. 

The LPS accepts scattered protons carrying a fraction of the beam momentum, xl, in the 
range xl ^ 0.4 and with < px ^ 1 GeV, where px is the transverse momentum of the proton 
with respect to the incoming beam direction. With the configuration installed in 1994 (S4, S5, 
S6), the resolution in xl is better than 0.4% at 820 GeV and the px resolution is about 5 MeV. 
The latter is less than the intrinsic transverse momentum spread in the proton beam at the 
interaction point (with rms of about 40 MeV horizontally and about 90 MeV vertically) due to 
the beam divergence of ~ 50 /irad in the horizontal and ~ 110 /zrad in the vertical plane. The 
LPS resolution is further discussed in section |3.1.2j . 



2.2.3 Reconstruction of an LPS track 

Tracks are reconstructed in stages, proceeding from individual hits to full tracks E2|. Noisy 



and malfunctioning channels are masked out and clusters of adjacent hit strips are searched 
for in each detector plane. Most clusters are one strip wide only (typically ~ 25% of the 
clusters have more than 1 strip). Track segments are then found independently in each detector 
assembly of six planes. As a first step, matching clusters in the two planes with the same strip 
orientation are combined. Candidate local track segments are then found by combining pairs 
of clusters belonging to different projections; when a pair of such clusters intersects within the 
region covered by the sensitive area of the detectors, a corresponding cluster in the remaining 
projection is searched for. In order to reduce the number of candidates, local track segments 
that traverse the overlap region of the detectors in the upper and the lower halves of the station 
are treated as one candidate. Finally, all hits belonging to a candidate (up to twelve for tracks 



4 



crossing the two halves, up to six otherwise) are used in a fit to find the transverse coordinates of 
the track at the value of Z corresponding to the centre of the station. The spatial resolution on 
these coordinates is about 30 /im. Figure f| shows the position of the reconstructed coordinates 
in the stations S4, S5 and S6 for a typical run. The regions with a high density of reconstructed 
hits in the overlap zone between the upper and the lower detectors correspond to tracks with xl 
close to unity. Lower xl tracks are deflected upwards and focussed horizontally onto a vertical 
line. For xl close to unity, this focus line is downstream of S6; it approaches S6 as xl decreases 
and reaches S6 for xl ~ 0.7. This explains the fact that for low xi tracks, the impact points in 
S5 and S6 tend to lie in a region which becomes narrower as the vertical coordinate increases. 

We distinguish two classes of events: those which are detected in all three of the stations 
and those which are detected in only two stations. In the latter, the interaction vertex position 
is required as a third point to measure the momentum. Tracks detected in three stations can 
be extrapolated backwards to Z = to also measure the transverse position of the interaction 
vertex. In both cases, coordinates reconstructed in pairs of different stations are first combined 
into track candidates and the track momentum is determined using the average ep interaction 
vertex with coordinates (X , Y ), found on a run-by-run basis with the sample of three-station 
tracks. Linear matrix equations relate the horizontal and vertical coordinates of the positions 
{hk,Vk) and slopes {h' k = dh k /dl, v' k = dvk/dl) of the track at each station to the position 
(Xo, Yq) and slope (X' , Yq) of the track at the interaction point. The coordinate along the 
beam trajectory is /. The positions (h k ,Vk) and slopes (h' k ,v' k ) are relative to the nominal 
beam position and direction at that value of I. The nominal beam crosses the interaction point 
(Z = 0) at X = Y = 0. For the horizontal direction one has: 

h k \_ m mi \fX \ fbo\ 
[h'J ~ { m 2 m 3 ) { X' ) [bj' 

An independent equation of the same form can be written for the vertical direction. The matrix 
elements are known functions of Xl which describe the beam optics including the effects 
of quadrupoles and drift lengths. The quantities (b ,bi), also functions of xl, describe the 
deflection induced by the dipoles and by the quadrupoles in which the beam is off axis; since 
the beam is taken as reference, they vanish as xl — > 1. 

Equation (|I|) and the corresponding one for the vertical direction can be written for a pair 
of stations (a, 6); upon eliminating the unknowns X' Q and Y Q ', one finds 



h b = M?(x L )h a + C£ b (x L ,X ), (2) 
v b = Mf(x L )v a + Cf(x L ,Y ), (3) 

where M ab and C ab are functions of the matrix elements and b{. These two equations 
are independent, apart from the common dependence on xl, and can be used to obtain two 
independent estimates of x^. If the values obtained are compatible, the pair of coordinates is 
retained CIS db candidate two-station track. 

As a final step of the pattern recognition, three-station track candidates are searched for 
using pairs of two-station candidates, e.g. one in S4, S5 and another in S5, S6. If a pair uses 
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the same hits in the station common to the two tracks, if the projections of the two tracks 
on the horizontal (non-bending) plane coincide and if the momenta assigned to each track are 
compatible, then the two candidates are merged in a three-station track candidate. Two-station 
and three-station candidates are then passed to a conventional track-fitting stage. 



A track \ 2 is defined as 



x 2 



E 



a 2 



(Xy - X ) 2 (Yy - Y ) 2 

n 2 n 2 ^ ' 

a X v a Y v 



where the sum runs over all clusters in all planes assigned to a track. Here Sj is the cluster 
position, Oi the uncertainty associated to it (which includes the effects of multiple scattering 
and the contribution of the cluster width; typical values range from 50 to 100 jum), (Xy,Yy) 
the interaction vertex coordinates in the X, Y plane, <Jx v and oy v the nominal widths of the 
vertex distribution; Si, a function of the five track parameters ip — (Xy,Yy,Xy,Y v ,Xi), is 
the predicted cluster position calculated from equation (|1]) and the corresponding one in the 
vertical direction; the quantities X v , Y v indicate the track slopes at the interaction vertex. The 
last two terms in eq. (^) constrain the track to the interaction vertex. This \ 2 is minimised 
with respect to the five track parameters, and the best track parameters, together with the 
error matrix, are determined. In the present analysis, for xl close to unity, the average value 
of x 2 / n df is ~ 1) where ndf is on average 7.3 for two-station tracks and 17.3 for three-station 
tracks. Three-station tracks are 60% of the total. 



2.2.4 Alignment of the LPS 

The alignment of the LPS relies on survey information for locating the detector planes in I 
and on high-energy proton tracks for locating them in h and v. The individual detector planes 
are first aligned within one station, then the relative alignment of the stations is determined, 
and finally the three stations S4, S5, S6 are aligned relative to the ZEUS detector. Typical 
accuracies in h and v are better than 20 /im. The actual path of the proton beam is also 
determined. These steps are described below. 

Tracks traversing the region in which the active areas of the detector planes in the upper and 
lower halves of a station overlap are used to align the detector planes within each half as well as 
to determine the position of the upper with respect to the lower half. With this procedure each 
plane is aligned independently; rotations of the detectors around the I axis are also determined. 

The relative alignment between the S4, S5, S6 stations in h is then found by exploiting the 
fact that tracks are straight lines in this projection. 

The only magnetic elements between S4 and S6 are the dipoles between S4 and S5 which 
deflect particles vertically. A sample of tracks with known deflection (i.e. known momen- 
tum) is thus necessary to align the stations relative to each other in v. This can be obtained 
independently of the LPS using the ZEUS calorimeter: 



x 



CAL 
L 



J2(E i+Pzi )/(2E p ), (5) 
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where the sum runs over all calorimeter cells, E is the energy measured in each cell and pz = 
E cos 6, with 9 the polar angle of each cell. The symbol E p denotes the incoming proton energy. 
Equation (Q) follows from energy and momentum conservation: + Pz)in = +Pz)out, 
where the sums run over the initial and final state particles, respectively, and J2(E + Pz)in = 
2E P . Events are selected with x^ AL > 0.99; these events have a clear peak in the xl spectrum 
as measured by the LPS, with very little background underneath. The relative positions of 
the stations are adjusted so that the peak appears at xl of unity. For the 1994 data, about 
20,000 events were used. The vertical alignment is finally checked by using events with elastic 
photoproduction of p° mesons - those discussed in the present paper - exploiting the fact that 
these events have scattered protons with xl very close to unity: xl can be written, for elastic 
p° photoproduction, as xl = 1 — (Q 2 + M 2 + \t\)/W 2 , where M p is the p° meson mass; for the 
sample used, the value of xl differs from unity by at most 0.2%. 

In order to align S4, S5 and S6 with respect to the proton beam line, tracks traversing 
all three stations are extrapolated to Z = 0, taking into account the fields of all traversed 
magnetic elements (mostly quadrupoles, as shown in Fig. p. The detectors are aligned with 
respect to the quadrupole axes by requiring that, independent of xl, the average position of the 
extrapolated vertex be the same as that measured by the central tracking detectors. At this 
point the detectors are aligned relative to the proton beam line and to the HERA quadrupole 
axes, and hence to the other components of ZEUS. About 40,000 three-station tracks were used 
for this procedure. 

Finally, the average angle of the proton beam with respect to the nominal beam direction 
is determined by using events of elastic photoproduction of p° mesons. For such events the 
transverse components of the scattered proton momentum balance on average those of the p° 
meson. The mean value of the sum of p x PS and p^™, and similarly for py, is set to zero 
by adding a constant offset to the fitted angle of the LPS tracks at the interaction vertex 
for all events. Here p x PS and Px TD indicate the X component of the proton momentum as 
measured by the LPS and of the p° momentum as measured by the CTD, respectively. This 
procedure defines the direction of the Z axis. Typical values of the beam offset are —15 /zrad and 
—100 /zrad in the horizontal and vertical directions, respectively, with respect to the nominal 
beam direction. The 1994 running period was split into three parts during which the beam tilt 
was relatively constant and the offset was determined for each part. Fig. |5| shows, separately 
for the X and Y projections, the sum of the proton and the p° transverse momenta after this 
correction, which is determined by requiring that both histograms be centred on zero. The 
width of the distributions is dominated by the intrinsic spread of the transverse momentum 
in the beam. The other (minor) contributions are the LPS and CTD resolutions and the fact 
that the transverse momentum of the scattered positron is not identically zero since Q 2 is not 
zero. Note that the effect of non-zero Q 2 is just to widen the distributions of Fig. |, not to 
shift them, since the X and Y components of the scattered positron momentum are centred 
on zero. In addition events with Q 2 > 0.01 GeV 2 contribute to the non-Gaussian tails. The 
Xl scale is not affected by this tilt correction. The sensitivity of the determination of t to the 
value of the tilt is weak, as discussed in section The effect of this correction is negligible for 
all quantities measured in the central ZEUS apparatus. 

As mentioned earlier, the detectors are in the retracted position between HERA runs; the 
positions of the pots (and hence of the detector planes) vary from one proton fill to the next 
by up to a few millimeters in Y (rarely in X) depending on the beam position and on the 
background conditions. Coordinate reconstruction can thus not be more accurate than the 
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reproducibility of the detector positioning system folded with the alignment accuracy This is 
monitored by the run-to-run dependence of the difference between the coordinates of the track 
impact point as measured by the detector planes in the upper and lower halves of a station for 
tracks in the overlap region. Note that this can be done since the alignment procedure described 
above is carried out using data from the whole running period, i.e. not on a run-by-run basis. 
The rms value of this difference is ~ 25 fim and is consistent with the specifications of the 
mechanics and commensurate with the detector resolution. 



3 Analysis 

3.1 Event selection 
3.1.1 Trigger 



ZEUS uses a three-level trigger system ]T2], For the present data, the trigger selected 

events with photoproduction of a vector meson decaying into two charged particles with no 
requirement that either the scattered positron or the scattered proton be detected. 

The first-level trigger required an energy deposit of at least 464 MeV in the electromagnetic 
section of RCAL (excluding the towers immediately around the beam pipe) and at least one 
track candidate in the CTD. Events with an energy deposit larger than 1250 MeV in the FCAL 
towers surrounding the beam pipe were rejected in order to suppress proton beam-gas events 
along with a large fraction of other photoproduction events. No requirements were made on 
the LPS information. 

At the second-level trigger, the background was reduced by using the measured times of the 
energy deposits and the summed energies from the calorimeter. 

The full event information was available at the third-level trigger; however, only a simplified 
reconstruction procedure was used. Tighter timing cuts as well as algorithms to remove cosmic 
muons were applied. One reconstructed vertex was demanded, with a Z coordinate within 
±66 cm of the nominal interaction point. Furthermore, the events were required to satisfy at 
least one of the following conditions: 

1. fewer than four reconstructed tracks and at least one pair with invariant mass less than 
1.5 GeV (assuming they are pions); 

2. fewer than six reconstructed tracks and no pair with invariant mass larger than 5 GeV 
(again assuming pions). 

Both sets of third-level triggers were prescaled by a factor six. Approximately 3 x 10 5 events 
were selected in this way, from an integrated luminosity of 898 ± 14 nb -1 (the luminosity 
corresponding to no prescale). 
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3.1.2 Offline requirements 



After performing the full reconstruction of the events, the following offline requirements were 
imposed to select elastic p° — > tt + tt~ candidates with a high momentum scattered proton: 

• Exactly two tracks in the CTD from particles of opposite charge, both associated with 
the reconstructed vertex. 

• The Z coordinate of the vertex within ±30 cm and the radial distance within 1.5 cm of 
the nominal interaction point. 

• In BCAL and RCAL, not more than 200 MeV in any EMC (HAC) calorimeter cell which 
is more than 30 cm (50 cm) away from the extrapolated impact position of either of the 
two tracks. This cut rejects events with additional particles, along with events with the 
scattered positron in RCAL. 

• One track in the LPS with 0.98 < xl < 1.02. This corresponds to a ±5cr window around 
xl = 1, for an xl resolution of 0.4%. As stated in section |2.2.4| , elastic photoproduction 
of p° mesons peaks at values of xl which differ from unity by less than 2 x 10~ 3 . This 
requirement is used to tag elastic events. 

• Protons whose reconstructed trajectories come closer than 0.5 mm to the wall of the 
beam pipe, at any point between the vertex and the last station hit, were rejected. This 
eliminates events where the proton could have hit the beam pipe wall and showered. In 
addition, it removes any sensitivity of the acceptance to possible misalignments of the 
HERA beam pipe elements. 



The value of the x 2 /ndf of the fit to the proton track (cf. section |2.2.3| ) less than 6 



The pion mass was assigned to each CTD track and the analysis was restricted to events 
reconstructed in the kinematic region defined by: 

0.55 < Mm < 1.2 GeV, 
0.27 < p T < 0.63 GeV, 
50 < W < 100 GeV, (6) 
Q 2 < 1 GeV 2 . 

The restricted range in the two-pion invariant mass M niT reduces the contamination from reac- 
tions involving other vector mesons, in particular from elastic (J) and ui production. The limits 
on pt, which is measured with the LPS, remove regions in which the acceptance of the LPS 
changes rapidly (cf . section |3.2|) . 



The photon-proton centre-of-mass energy W and the mass were determined from the 
momenta of the two pions 0. Energy and momentum conservation relate the photon energy, 
Ej, to the two-pion system energy E n7T and longitudinal momentum pzmr by 2E 7 « (E wn — 
Pztttt), under the assumption that the positron emits the virtual photon with zero transverse 
momentum. Therefore W 2 ~ AE 1 E p « 2(E nn —pzmr)E p . From the Monte Carlo study discussed 
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in the next section, the resolution on W has been found to be about 2 GeV; that on M 7T7T is 
about 30 MeV. 

The combination of the trigger requirements and the offline cuts, excluding that on Q 2 , 
limits Q 2 to be less than »s 4 GeV 2 . However, unlike previous ZEUS analyses of untagged 
photoproduction events, for the present data Q 2 was determined event by event. By exploiting 
the transverse momentum balance of the scattered positron, the tt + it~ pair and the scattered 
proton, one obtains pj>, the transverse momentum of the scattered positron. The variable 
Q 2 was then calculated from Q 2 = (p^) 2 /(l — y), where y is the fraction of the positron 
energy transferred by the photon to the hadronic final state, in the proton rest frame; it 
was evaluated as y ~ W 2 /s, where y/s is the ep centre-of-mass energy. Figure || shows a 
scatter plot of the reconstructed and the generated value of Q 2 for the sample of Monte Carlo 
events used to evaluate the acceptance (cf. section |3~2]) ; the line shows the expected average 
relationship between these two quantities assuming a beam transverse momentum distribution 
with <j Vx = 40 MeV and a PY = 90 MeV. At small values of Q 2 , the resolution on Q 2 is dominated 
by the beam transverse momentum spread; it is about 100% at 3 x 10~ 2 GeV 2 , 40% at 0.1 GeV 2 
and 20% at 1 GeV 2 . 

The final sample contains 1653 events. Figure |7| shows the M n7T , W, px, Py, Pt and xl dis- 
tributions after the offline selections; the variables px and py denote the transverse components 
of the outgoing proton momentum with respect to the incoming beam axis and p\ = p\ + Py- 
The invariant mass plot is dominated by the p° peak. The shape of the px spectrum, with two 
well separated peaks, is a consequence of the fact, discussed earlier, that events with xl close 
to unity only populate a narrow region of the detectors at v ~ 0. As discussed earlier, elastic 
p° photoproduction peaks at xl — 1 to within < 2 x 10 -3 ; the width of the xl distribution 
in Fig. [7]f shows that the resolution of the LPS in xl is ~ 0.4%. For the same events, Fig. |] 
shows the scatter plot of the reconstructed X and Y components of the proton momentum. 
For the present measurement, only the pt region between the dashed vertical lines in Fig. ^e 
was used. The variable t = (p — p') 2 , where p and p' are the incoming and the scattered proton 
four-momenta, respectively, can be evaluated as follows: 

t= (P - P'f « [l + (M 2 p /p 2 T )(x L - l) 2 ] , (7) 

where M p is the proton mass and terms of order (M 2 /p' ) 2 or higher are neglected. For the 
present events, which have xl ~ 1, the approximation t ~ —Pt/ x l ~ —p\ was used. 

Finally, Fig. |9| shows the distribution of the reconstructed values of Q 2 . As discussed above, 
at small values of Q 2 the intrinsic spread of the beam transverse momentum dominates. The re- 
quirement that Q 2 be less than 1 GeV 2 removes 7 events. The median Q 2 of the data, estimated 
with the Monte Carlo study discussed in the next section, is approximately 10~ 4 GeV 2 . 



3.2 Monte Carlo generators and acceptance calculation 



The reaction ep ep°p was modelled using the DIPSI generator, which was shown to 
reproduce the ZEUS p° photoproduction data 0. The effective W dependence of the jp cross 
section for the events generated was of the type a oc W ' 2 . The t distribution was approximately 
exponential with a slope parameter of 9.5 GeV -2 . The two-pion invariant mass, M^ n , was 
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generated so as to reproduce, after reconstruction, the measured distribution. The angular 
distribution of the decay pions was assumed to be that expected on the basis of s-channel 
helicity conservation |24| . 



The simulated events were passed through the same reconstruction and analysis programs 
as the data. In Figures [7] and ^| the distributions of the reconstructed data (not corrected 
for acceptance) over M^, W, px, Py, Pt, xl and Q 2 are compared with those obtained for 
the reconstructed Monte Carlo events. The Monte Carlo is in reasonable agreement with the 
data. Figure |TU| a shows the overall accept ctncG £is db function of t, obtained using DIPSI. The 
acceptance includes the effects of the geometric acceptance of the apparatus, its efficiency and 
resolution and the trigger and reconstruction efficiencies. Since the detector planes cannot be 
positioned in the beam, the acceptance vanishes at small values of t. Conversely, in the PxiPy 



region covered by the detectors, the acceptance of the LPS is large, as shown in Fig. |10|b, which 
shows the geometric acceptance of the LPS alone, irrespective of the acceptance of the rest of 
the ZEUS apparatus. The region of LPS geometric acceptance larger than 95% for both px > 



and px < maps into that of 0.25 <Pt^j 0.65 GeV; as discussed in section |3.1.2j , events outside 



this region are not used in the present analysis. For events with elastic p° photoproduction, 
the geometric acceptance of the LPS, averaged over azimuth, is approximately 6%. 

As discussed in the next section, in order to estimate the contamination from the reaction 
ep — > ep°Xx, where is a hadronic state of mass resulting from the diffractive dis- 
sociation of the proton, the PYTHIA generator [2£| was used. A cross section of the form 



d 2 a/dtdM 2 N oc e-W/Mfr, with b = Q GeV~ 2 , was assumed; the value of Mjy ranged between 
M p + M n (where M p is the proton and M n the pion mass) and a maximum fixed by the condition 
Mfj/W 2 < 0.1 [p^j . The p° decay angular distributions were assumed to be the same as those 
of the elastic events. 



3.3 Backgrounds 



After applying the selection criteria described in section 3.1, the data contain small contribu- 



tions from various background processes to the reaction ep — > eiT + 7i p: 

• Beam-halo tracks observed in the LPS may overlap with events in which a p° is seen by 
the ZEUS central apparatus. The term beam-halo refers to protons with energy close to 
that of the beam originating from interactions of beam protons with the residual gas in 
the pipe or with the beam collimators. Such tracks are completely uncorrelated with the 
activity in the central ZEUS apparatus; therefore, any sample of events selected without 
using the LPS information contains the same fraction, ehaio, of random overlaps from halo 
tracks within the LPS acceptance. This fraction was found to be ehaio = 0.25 ± 0.03% 
by analysing events of the type ep — > eX at Q 2 > 4 GeV 2 in which the virtual photon 
diffractively dissociates into the state X. For these events one can measure X and the 
scattered positron in the calorimeter; in addition a proton track is looked for in the LPS. If 
one is found, the event is fully contained and its kinematics is thus overconstrained: most 
beam-halo events appear to violate energy-momentum conservation and can therefore be 
identified. 

The contamination of the present sample (after the requirement of a good LPS track) 
can be obtained as (eh a ioN no LPs) / Nlps = 5.0% ± 0.6% (stat.), where N no LPs indicates 
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the number of events found by applying all cuts except for the requirement that a track 
be detected in the LPS, and N LPS = 1653 is the number of events after all cuts. These 
events were not removed from the present sample, but their effect on the measurement of 
the t slope is small, as discussed in section §. 

In the reaction ep —>■ ep°Xjv, the proton diffractively dissociates into a hadronic state 
Xn which may escape detection by the central detector. The debris of may con- 
tain a proton within the LPS acceptance and with 0.98 < xl < 1-02: such events are 
indistinguishable from elastic p° production. 

In order to evaluate the contamination from such events, the cut on xl was removed; 



Fig. [11] shows the Xl spectrum thus obtained, not corrected for acceptance. The sum of the 
reconstructed xl distributions from DIPSI and PYTHIA was fitted to this spectrum with 
the normalisations of the simulated distributions as free parameters of the fit. The fit gives 
an acceptable description of the data, as shown in Fig. [TTJ. The resulting contamination 
of proton-dissociative events for xl > 0.98 is 0.21% ±0.15% (stat.), a major improvement 
with respect to 11% ± 1% (stat.) ± 6% (syst.) in the earlier ZEUS result || which did not 
use the LPS. 

The contaminations from elastic production of uo and <fi mesons were estimated in [0 to 
be (1.3 ± 0.2)% and (0.3 ± 0.1)%, respectively. 

Contamination from positron beam-gas and proton beam-gas events was studied by using 
the unpaired bunches event samples to which all the cuts described above were applied. 
No event passed the cuts, indicating a negligible contamination. 



4 Results 

The differential cross section da/dM^ for the process 'jp — > 7r + 7r~p was evaluated in the 
kinematic range 0.55 < M m < 1.2 GeV, 50 < W < 100 GeV, Q 2 < 1 GeV 2 and 0.073 < \t\ < 
0.40 GeV 2 . In each bin the cross section was obtained as 



^ Cy "' (8) 



where N 7T + 7T - is the number of observed events in the bin, L is the integrated luminosity, a is 
the overall acceptance in the bin, and $ = 0.0574 is the photon flux factor, i.e. the integral of 
equation (5) in ref. || over the measured W and Q 2 ranges of this measurement. The factor 
Chaio = 0.950 ± 0.006 (stat.) corrects for the beam-halo contamination discussed in section 



The effects of positron initial and final state radiation and that of vacuum polarisation loops 
were neglected; the effects on the integrated cross section were estimated to be smaller than 
The effects on the shape of the M T7r and t distributions are expected to be negligible. 



The small residual contaminations from proton dissociative p° production and elastic uo and 
production, discussed in the previous section, were not corrected for. 

Figure |T2"| shows the differential cross section da/dM n7T in the interval 0.55 < M n7T < 1.2 GeV, 
0.073 < \t\ < 0.40 GeV 2 for (W) = 73 GeV. The mass spectrum is skewed, as previously 
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observed both at low energy and at HERA. This can be understood in terms of resonant and 
non- resonant 7r + 7r~ production ||29|| , and their interference J30|. The spectrum was fitted using 



expression (11) of ref. 0]. The results for the total, resonant and interference terms, as obtained 
in the fit, are indicated in the figure. The fraction of the resonant to the total contribution in 
the measured range was found to be c res = 0.91 ± 0.04 (syst.). The uncertainty was evaluated 
by repeating the fit with the various functional forms discussed in 0. In the contribution 
of the resonant term was found to vary from 86% for \t\ = 0.01 GeV 2 to 95% for \t\ = 0.5 
GeV 2 . No t dependence of c res was assumed here, except in the evaluation of the systematic 
uncertainty (see below). 

The differential cross section da/dt for the reaction 7^ — > p°p was obtained similarly to 
da/dM n7T , but in addition the correction factor c res , just discussed, was applied. Figure [i~3| shows 
the result in the interval 0.073 < \t\ < 0.40 GeV 2 , 0.55 < M ww < 1.2 GeV for (W) = 73 GeV. 
The data were fitted with the function 



dt~ A 6 ' 



(9) 



the result of the fit is shown as a straight line on Fig. [13|. The fitted value of the slope parameter 
b is 



b = 9.8 ± 0.8 (stat.) ± 1.1 (syst.) GeV" 2 . (10) 

The result is consistent with b = 9.9±1.2 (stat.)±1.4 (syst.) GeV~ 2 obtained in 0J for the range 
60 < W < 80 GeV, Q 2 < 4 GeV 2 and |t| < 0.5 GeV 2 using a fit of the type Aexp (-b\t\ + ct 2 ). 
For both the present data and those of ref. 0J, (W) ~ 70 GeV. 

The measured differential cross section was integrated over the range 0.073 < \t\ < 0.40 GeV -2 , 
yielding o = 5.8 ± 0.3 (stat.) ± 0.7 (syst.) /ib, again at (W) = 73 GeV and for 0.55 < M wn < 
1.2 GeV. The result can be extrapolated to the mass range 2M n < M nn < M p + 5r , as 
in 0, using the fit to the mass spectrum described earlier (here Tq is the p° width); this yields 
o = 6.3 ± 0.3 (stat.) ± 0.8 (syst.) /ib, where no uncertainty was assigned to the extrapolation. 
If our previous result M is integrated in the t range covered by the present data, using the 
published results of the fit with the function Aexp (— b\t\ + ct 2 ) (table 5 of ref. @, left column), 
one finds a = 6.7 ± 1.1 (syst.) /ib, in good agreement with the present result; only the system- 
atic uncertainty is given since it is dominant. This uncertainty was obtained by scaling the one 
published in M for the cross section measured over the range \t\ < 0.5 GeV 2 by the ratio of the 
cross sections for the present t range and for \t\ < 0.5 GeV 2 . 

The major sources of systematic uncertainty on b and a are the acceptance determination 
and the background contamination, the former being dominant. Table [I] lists the individual 
contributions. In the following we discuss them in detail. 



1. In order to estimate the uncertainty due to the acceptance, the analysis was repeated 
varying the requirements and procedures as listed below. 
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Contribution 


|A6/6| 


|Acr/cr| 


Integrated luminosity 


- 


1.5% 


Acceptance: trigger efficiency 


- 


9% 


Acceptance for pion tracks 


< 1% 


1% 


Acceptance for proton track 


7% 


6% 


Acceptance: sensitivity to binning in t 


2% 


- 


Acceptance: unfolding of beam transverse momentum spread 


7% 


- 


Acceptance: sensitivity of p beam angle 


3% 


1% 


Background: beam-halo 


4% 




Procedure to extract the resonant part of the cross section 


1.6% 


4% 


Background due to elastic uo and <fi production 




1% 


Radiative corrections 




4% 


Total 


11% 


12% 



Table 1: Contributions to the systematic uncertainty on b and a. 

(a) For the pion tracks in the central detector: 

• The pseudorapidity rj = — In tan (0/2) of each of the two tracks was restricted 
to the range \rj\ < 1.8, thereby using only tracks which have traversed at least 
three superlayers in the CTD. 

• The radial distance of the vertex from the beam axis was required to be less 
than 1 cm. 

In both cases the changes are small; by summing them in quadrature one finds 
\Ab/b\ = 0.2% and \Aa/a\ = 1%. 

(b) For the proton track in the LPS: 

• The maximum allowed value of y 2 /ndf for the reconstructed proton track was 
reduced from 6 to 2. 

• The minimum distance of approach of the proton trajectory to the beam pipe 
was increased from 0.5 mm to 1.5 mm. 

• Events with p x > and with p x < were analysed separately, as a check of 
possible relative rotations of the stations. 

• The data were divided into a "large acceptance" and a "low acceptance" sample 
depending on the position of the LPS stations, which, as discussed above, varied 
slightly from run to run. 

By summing the individual contributions to Ab/b in quadrature, independently of 
their sign, \ Ab/b\ = 7% is obtained. The corresponding uncertainty on a is |A<j/a| = 
6%. 

(c) The sensitivity of the result on b to the binning in t was studied by reducing bin 
sizes by up to 20%; the bin edges were moved by up to one fourth of the bin size. 
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The largest effect was 2% for |A6/fe|. 

(d) As discussed earlier, t has been obtained as — p 2 -., with p\ = p\ + p\ the trans- 
verse momentum of the scattered proton with respect to the incoming beam axis. 
Since the incoming proton beam has an intrinsic transverse momentum spread of 
a Px « 40 MeV and a PY ~ 90 MeV, which is much larger than the LPS resolution 
in transverse momentum, the measured value of t is smeared with respect to the 
true t. The Monte Carlo simulation takes into account the proton beam transverse 
momentum spread. The acceptance corrected t distribution is thus corrected also 
for this effect. 

The following alternative approach to account for the effect of the transverse momen- 
tum spread of the beam has also been followed. Assuming that the true t distribution 
has the form given by equation (||), the measured p\ distribution can be expressed 
as a convolution of equation (^) and a two-dimensional Gaussian distribution repre- 
senting the beam transverse momentum distribution, with standard deviations o Px 
and a PY . Unfolding the contribution of the beam transverse momentum spread from 
da /dp 2 -, provides an alternative evaluation of da/dt. In this case one first measures 
the distribution of p\ without making any correction for the effects of the beam 
intrinsic spread, thereby exploiting the good resolution of the LPS on the trans- 
verse momentum. In a second stage, the effect of the beam spread is unfolded. If 
a px = 40 MeV and a PY = 90 MeV (as seen in the data, cf. Fig. |), then |A6/6| = 7%, 
with only a weak dependence on the values of a Px and a PY . 

(e) The sensitivity to the determination of the proton beam angle (cf. section |2.2.4j ) 
was evaluated by systematically shifting p? by 10 MeV. This amount is twice the pt 
resolution of the LPS and corresponds to > 5 times the uncertainty on the means of 
the distributions of Fig. |5|. The corresponding variations of b and a were | Ab/b\ = 3% 
and \Aa/a\ = 1%. 

The differences between the values of b obtained in cases (a) to (e) and that obtained 
with the standard analysis were summed in quadrature, yielding \Ab/b\ = 10.5% and 
\Aa/a\ = 6%. 

2. Effect of background contamination. 

(a) As mentioned above, no correction was applied for a possible t dependence of the 
background. The only significant background is the halo. If the assumption is made 
that the halo contribution (5.0%±0.6%) has a distribution of the type exp (— bh a io\t\), 
then |Afe/6| < 4% when bhaio is varied between 5 and 15 GeV -2 ; this range of 
variation is consistent with estimates of bhaio based on the ep — > eXp events at 
Q 2 > 4 GeV 2 discussed in section [O. 



(b) If the t dependence of c res evaluated in is assumed for the present data, the slope 
changes by Ab/b = —1.6%. 



The latter two contributions were also added quadratically to the systematic uncertainty, yield- 
ing a total systematic uncertainty of 11% on b, dominated by the LPS acceptance and the effect 
of the beam transverse momentum spread. 

The total systematic uncertainty on a is 12%, which includes, in addition to the contributions 
detailed above, the uncertainty on the luminosity (1.5%), that on the trigger efficiency (9%) 
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and that related to the extraction of the resonant part of the cross section. The estimated 
background due to elastic u and production, as well as the upper limit of the correction 
for radiative effects have also been included. The systematic uncertainty on a is dominated 
by contributions not related to the LPS (11%); the uncertainty on the LPS acceptance is 6%, 
which has only a small effect on the total uncertainty when summed in quadrature with the 
other contributions. 



5 Conclusions 

The Leading Proton Spectrometer of ZEUS is a large scale system of silicon micro-strip detectors 
which have been successfully operated close to the HERA proton beam (typically a few mm) 
by means of the "Roman pot" technique. It measures precisely the momentum of high energy 
scattered protons, with accuracies of 0.4% for the longitudinal and 5 MeV for the transverse 
momentum. 

As a first application, the cross section, the M niT and the t dependences of the reaction 
IP ~* P°P have been measured in the kinematic range Q 2 < 1 GeV 2 , 50 < W < 100 GeV, 
0.55 < < 1.2 GeV and 0.073 < \t\ < 0.40 GeV 2 . Elastic events were tagged by demanding 
that xl he larger than 0.98, i.e. that the scattered proton carry at least 98% of the incoming 
proton momentum. For the first time at these energies, t was measured directly. Compared 
to our previous analysis, the present technique based on the use of the LPS eliminates the 
contamination from events with diffractive dissociation of the proton into low mass states. 

In the range 0.073 < \t\ < 0.40 GeV 2 , the differential cross section da/dt is described by 
an exponential distribution with a slope parameter b = 9.8 ± 0.8 (stat.) ± 1.1 (syst.) GeV -2 . 
The systematic uncertainty is dominated by the uncertainty on the LPS acceptance and the 
effect of the intrinsic transverse momentum spread of the beam. In the measured t and M n7T 
intervals, the integrated p° photoproduction cross section at (W) = 73 GeV was found to be 
5.8 ± 0.3 (stat.) ± 0.7 (syst.) fib, consistent with our previous measurement obtained in a 
slightly different kinematic range. 
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Figure 1: Side view of the outgoing proton beam line for Z between 5 and 50 m (bottom) 
and between 50 and 90 m (top). The small drawing below the Z axis is to scale; in the 
large drawing different scales are used for the longitudinal and transverse directions. Magnetic 
elements labelled BH, BS and BT are horizontally bending dipoles; vertically bending dipoles 
are indicated as BZ, BY and BU. Quadrupoles are labelled as QL, QK, QC, QB, QS, QR; the 
magnets upstream of the septum magnet BS are common to the proton and positron beam 
line. The positions of the LPS stations SI through S6 are also shown. The centre of the ZEUS 
detector is at Z = 0. 
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Figure 2: Schematic layout of a station (like S4, S5 or S6). A) During beam filling and ramping, 
the detector planes (labelled "Silicon detector" ) are kept outside of the pots and the pots are 
placed far from the beam. The zig-zag lines indicate the bellows. B) The detector planes are 
inside the pots and the pots are being moved towards the beam. Note the elliptical profile of 
the fronts of the pots (X-Y view), which matches the cutout of the detector planes. C) When 
taking data, the pots are fully inserted and the detector planes in the upper and lower half of 
the station partially overlap in the transverse plane. 
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Figure 3: Simplified diagram of the HERA proton beam optics relevant for the LPS. The 
top view (labelled "horizontal") shows the magnetic elements relevant to the SI, S2 and S3 
spectrometer as thick lines. Note the half quadrupole between SI and S2. The side view 
(labelled "vertical") shows the magnetic elements relevant to the S4, S5 and S6 spectrometer 
as thick lines. 
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Figure 4: Positions of the reconstructed track impact points in S4, S5 and S6. For each plot 
the origin of the reference frame coincides with the position of the nominal proton beam at 
the value of Z corresponding to the centre of the station. The continuous lines approximately 
indicate the sensitive region of the detector planes. 
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Figure 5: Sum of the X (a) and Y (b) components of the proton momentum as measured by the 
LPS and of the p° momentum as measured by the CTD. The histogram was obtained with the 
Monte Carlo simulation discussed in section [O . The continuous line is the result of a Gaussian 
fit to the data. The fitted values of the standard deviations are 45 and 102 MeV, respectively 
They are dominated by the spread of the transverse momentum in the beam. The other minor 
contributions are the LPS and CTD resolutions and the fact that Q 2 is not exactly zero; events 
with Q 2 > 0.01 GeV 2 contribute to the non-Gaussian tails. 
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Figure 6: Scatter plot of reconstructed versus generated values of Q 2 for Monte Carlo events. 
The continuous line shows the expected average relation assuming a spread of the beam trans- 
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Figure 7: Observed distributions for (a) M T7T , (b) W, (c) p x , (d) py, (e) p T and (f) of the 
reconstructed data (points) and the reconstructed Monte Carlo events (histogram). The distri- 
butions are not corrected for acceptance. The vertical bars indicate the statistical uncertainties. 
The dashed lines in (e) show the limits of the pt region used. 
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Figure 8: Scatter plot of the X and Y components of the scattered proton momentum as mea- 
sured by the LPS for the accepted events. The continuous curves correspond to the indicated 
values of \t\ (in GeV 2 ), which were used as limits of the bins in Fig. [13|. 
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Figure 9: Q 2 distributions of the reconstructed data (points) and the reconstructed Monte 
Carlo events (histogram). The distributions are not corrected for acceptance. Only the region 
Q 2 < 0.1 is shown. The vertical bars indicate the size of the statistical uncertainties. 
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Figure 10: (a) Acceptance as a function of £; it includes the effects of the geometric acceptance 
of the apparatus (ZEUS and LPS), of its efficiency and resolution and of the trigger and 
reconstruction efficiencies. The vertical bars indicate the size of the statistical uncertainties, 
(b) Purely geometric acceptance e for the LPS alone for xl = 1 tracks as a function of px and 
Py- 
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Figure 11: Reconstructed spectrum for the data (points). All selection cuts were applied 
except that on x^. The distribution is not corrected for acceptance. The sum of the recon- 
structed Xi distributions produced with the DIPSI and PYTHIA generators was fitted to this 
spectrum with the normalisations as free parameters of the fit. The result of the fit is shown 
as a histogram. The DIPSI generator simulates the elastic reaction ep — > ep°p, and PYTHIA 
the proton dissociative reaction ep — > ep^X^- The contribution of PYTHIA is shown as the 
hatched area. The vertical bars indicate the size of the statistical uncertainties. 
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Figure 12: Differential cross section da/dM nn for the reaction 7£> — > 7r + 7r~p, in the region 
0.55 < M,,-^ < 1.2 GeV and for (W) = 73 GeV. The vertical bars indicate the size of the 
statistical uncertainties only. The lines indicate the result of the fit using expression (11) of 
ref. §• 
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Figure 13: Differential cross section da/d\t\ for elastic p° photoproduction, 7p — > p°p, at 
(W) = 73 GeV in the region 0.073 < \t\ < 0.40 GeV 2 , 0.55 < M m < 1.2 GeV. The vertical bars 
indicate the size of the statistical uncertainties only. The line is the result of the fit described 
in the text. 
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